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ABSTRACT 


This thesis examines the use of heat for producing clinical lesions 
in tissue and presents the design and analysis of a resistively heated 
Surgical probe. The probe temperature is accurately maintained and 
measured by using a Wheatstone bridge. The probe was embedded in an 
agar-water test medium and the temperature field emanating from the 
probe was measured using liquid crystal tapes. The experimental results 
compare within approximately 10% of a two-dimensional numerical solution. 
A one-dimensional theoretical model is developed and results are 
presented to show the relative effects of metabolism and blood flow on 


the temperature field. 





el. 


LF 


IV. 


wel. 


Weel; 


TABLE OF CONTENTS 


Mine ND UGIR = eeeenee aoe be tte. -- =e 
SUMMAR VAC are MING EGHIEOUES) -----=-----=------------— 
PEE DUR EGIMGURR ENT (de) METHOD: =---=-2=--=2-=-=-~=----- 
B. RADIO-FREQUENCY (rf) METHOD ----------------------- 
STIG RESISTANCE SHEAING: METHOD) =-------2==----=--.- 
MGM SONS ANCE: HEAT INGEPROBE (== --—-=222-=-==—2---22--- 
DME EQ RMS Sos ae = Sone tae asa ne see ss caee ase. 
PRE ESIC eee = 2 ens ene ase ee een ee ones 
PRIENESE AUR ERGONTIRULNs == seen === eee 8 
BORNE APPARATUS 222---2--22=- 22-2 -5--- 22 2c--—- 
A, TES PRED es eee 


B. LIQUID CRYSTALS AND TEMPERATURE MEASUREMENT ------- 


C. TEST CELLS ---------------------------------------- 
EXPERIMENTAL PROCEDURE -------------------------------- 
THEORETICAL ANALYSIS ---------------------------------- 
A. ONE-DIMENSIONAL MODEL ----------------------------- 

1. Metabolic and Blood Flow Effects Negtected ---- 

2. Metabolic and Blood Flow Effects Included ----- 
B. TWO-DIMENSIONAL MODEL ----------------------------- 
es 
A. EXPERIMENTAL RESULTS ------------------------------ 
B. ONE-DIMENSIONAL MODEL ----------------------------- 


1. Effect of Blood Flow on the Temperature 
FIC@Id ------ nnn nnn nn ne nnn nnn nn ene nen n ene 


2. Effect of Increased Time on the Temperature 
Field ----------------------------------------- 


1] 
13 
As 
14 
iis 
18 
18 
ng 
2 | 
ZS, 
Za 
30 
33 
ey 
4] 
4] 
4] 
43 
45 
47 
47 
62 


62 


62 





; 
‘ 


Saas 


Pe ORE owen AL Obie 9s —————~—=-—==—— == 66 


VIII. CONCLUSIONS AND RECOMMENDATIONS ----------------------- 68 
PMEGONGIISIONS i222" =2-=-2----=--~-~-~ SMa... aM 68 
ERAGON END ONSW= == eames a2 2a 5----------.--.-- 69 

APPENDIX A: ONE-DIMENSIONAL MODEL SOLUTION ----------------- 70 

APPENDIX B: ERROR ANALYSIS OF EXPERIMENTAL RESULTS --------- 74 

Ber ERED ROGRMI Milli == 252 sse-22252=--252-2--5__ 2225.44... 94 

Pe IER MER OGRE aie se= == -- eee oe 2s __--..- 105 

CANE REROGRAN| <9 “e2s22--2-2ss-es—5->-----2-5e2-5---- 2s 108 

MERRUTER PROGRAM 4 --=----------------2-.---_----=---..---- 111 

UCRUNERMeROGRAM 5) 2 =-2=2=---- 2222 a2 ee Tee ___- 114 

OMe BIEREDROGRIM a6. 2o=~o2----22-— 8-8 8-85 =e 117 

2S TOC ee 120 

Meee MRDISTRIGUTION LIST --=--2--222-----=---------s----=--- 122 

er UI Cvicis = o-  en e Se SSS ee... __-.- 123 


se 





Ii. 


Welel 


Iv 


VI. 


Wer. 


VIII. 


I. 


LIST OF TABLES 


Resistivity and the temperature coefficient of 


resistivity for various metals ----------------------- 


Theoretical values for 6 for various values of 
nondimensional location (R) and values of t 
between 0.1 and 0.9: one-dimensional solution, 


ee ee 


Theoretical values for 6 for various values of 
nondimensional locations (R) and values of rt 
between 1.0 and 9.0: one-dimensional solution, 


0 ee ee ee 


Theoretical values for 6 for various values of 
nondimensional location (R) and values of t 
between 10.0 and 100.0: one-dimensional solution, 


ee eee ee 


Theoretical values for 6 for various values of 
nondimensional locations (R) and values of t 
between 0.1 and 0.9: one-dimensional solution, 


PO oe ae eee ea eee een e ede 


' Theoretical Values for 6 for various values of 


nondimensional location (R) and values of rt 
between 1.0 and 9.0: one-dimensional solution, 


Pe Oe ene een seen a ene meee eae Sere 


Theoretical values for 6 for various values of 
nondimensional locations (R) and values of t 
between 10.0 and 100.0: one-dimensional solution, 


POM ee ee oer e nese eee 


Theoretical values for 6 for various values of 
nondimensional location (R) and values of rt 
between 0.1 and 0.9: one-dimensional solution, 


Pe OS tsae aoe a ae ws eee eee ee eee 


Theoretical values for 6 for various values of 
nondimensional locations (R) and values of rt 
between 1.0 and 9.0: one-dimensional solution, 


§ S005 s secs assesses sa ses se se aes 


Theoretical values for 6 for various values of 
nondimensional location (R) and values of t 
between 0.1] and 0.9: one-dimensional solution, 


3 eee ee ee nee eee 


20 


76 


77 


78 


i 


80 


8] 


82 


83 


84 





XI. Theoretical values for 6 for yarious values of 
nondimensional locations (R) and values of rt 
between 1.0 and 9.0: one-dimensional solution, 
8 = 1.0 ---------------------------------------------- 


XII. Theoretical values of 6 for various values of 
nondimensional time (t) and nondimensional Tocation 
(R): two-dimensional solution, 8 = 0.0, X=I1R --------- 


XIII. Theoretical values of 6 for various values of 
nondimensional time (+t) and nondimensional location 
(R): two-dimensional solution, 8 = 0.0, X=5R --------- 


XIV. Theoretical values of 6 for various values of 
nondimensional time (+) and nondimensional location 
(R): two-dimensional solution, 8 = 0.0, X=8.5 -------- 


XV. Theoretical values of 6 for various values of 
nondimensional time (rt) and nondimensional location 
(R): two-dimensional solution, g = 0.0, X=11.5 ------- 


XVI. Theoretical values of 6 for various values of 
nondimensional time (t) and nondimensional location 
(R): two-dimensional solution, 8 = 0.0, X=13.5R ------ 


XVII. Experimental data ------------------------------------ 





al: 
ee 
ee 
14. 


OF 
We. 
Ise 
no. 
Zi. 
Zl 
le 
“a3 
24. 


Construction of probe 


Completed probe 


Wheatstone bridge 


LIST OF ILLUSTRATIONS 


Temperature control circuit 


Probe calibration curve 


Test cells 


Removable support rings for liquid crystal tape ------ 


Top view of circular test cell 


Overall view of apparatus 


Mathematical model of probe 


TRUMP nodal network 


Radial 
Radial 
Radial 
Radial 
Radial 


temperature 
temperature 
temperature 
temperature 


temperature 


field for 
field for 
fel die kor, 
field for 
field for 


Axial temperature field for 


Axial 
Axial 
Axial 


Axial 


Rvs. 


c 


R_ vs. 


C 


Axial position (X) vs. 
of t and ae 0.285 


temperature 
temperature 
temperature 
temperature 
+ for Oo = 


T For. 6. = 


field for 
field for 
field for 
field for 
0.200 and 
0-235 and 


C 


B 


B 


R for various values 


gae gae@taeete@tFteee Sez se eee SS SPP ae eae = se ae SP ae ee te FF ae a= = 


22 
23 
24 
26 
27 
34 
35 
38 
40 
42 
46 
4g 
49 
50 
51 
52 
53 
54 
55 
56 
57 
59 
60 


6] 





Ere Pere isan 
| Peete eel 


Bos 
26. 
2/. 
we. 


6vs. R fort 


=—or) ana 6 = 0.0, 0.1, O:Smand 1.0 ---- 


6 vs. t for various values of R and B = 0.0 ---------- 


6 vs. t for various values of 8 and R = 2.5 --------- - 


Axial position (X) vs. R for various values 


of t and 6 


0.300 


64 
65 


67 





NOTATION 


Description 


Tissue thermal diffusivity 
Temperature coefficient of resistance 
Tissue specific heat 

Blood specific heat 
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Heating rate 

Radial location in tissue 
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Matabolic heat generation rate 
Time 

Tissue temperature 

Arterial temperature 
Critical temperature 
Initial: tissue temperature 
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I. INTRODUCTION 


A variety of surgical techniques are presently used to produce 
clinical lesions. In addition to the familiar mechanical or "scalpel" 
method, surgical procedures employing chemical, ultrasonic, focused 
X-rays, radiation, and cooling and heating agencies are used to destroy 
descrete regions of tissue, normally nervous tissue. The choice of 
methods depends largely on the nature of the clinical problem. Ideally, 
the surgeon desires a technique that possesses the following character- 
emies: 1) Safety, 2) Reversibility, 3) Reproducibility, 4) Sharp 
delimitation, 5) Hemostasis, 6) Flexibility, 7) Simplicity, and 8) 

Speed of application. At times, it 1s not possible to meet all of the 
above surgical ideals. In such cases, the choice of procedures 1S 
based largely on an optimization of the ideals, with safety being of 
paramount importance. Gengler [10] presents a summary of the various 
methods, complete with the advantages and disadvantages of each. 

In this thesis, the use of heat for producing lesions is investigated. 
In particular, the design of a resistively heated surgical electrode, 
henceforth referred to as a probe, is presented. Tissue destruction is 
accomplished by using the probe as a heating element to raise the 
temperature of a predetermined region of tissue above a critical level, 
normally taken to be 55°C [2]. 

The unique feature of this probe is that its surface temperature 
can be accurately measured and controlled by using a Wheatstone bridge. 
As a result, the temperature field in the medium surrounding the probe 
can be analytically predicted by solving the heat equation prior to 


applying the probe. This is quite significant when one considers that 


1] 
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the volume of tissue destroyed depends solely on the local tissue 
temperature. 

Temperature fields determined experimentally with the probe 
embedded in a 0.30% agar - 99.7% water test medium compare within an 
accuracy of 13% of those predicted using the analytical model. A novel, 
relatively new, temperature sensing device was used in the experimental 
phase of the study. Liquid Crystals [12], a material that changes color 
over a known and well defined temperature range, were employed to 
obtain a two-dimensional, visual display of the transient temperature 
field which developed around the probe. By using two separate orien- 
tations of the liquid crystals, a three-dimensional picture of the 
temperature field was inferred. In addition to the highly desirable 
visual aspect of the liquid crystal material, it also produces a minimum 
disturbance in the test medium. 

The effects of blood flow and metabolism on the size of the lesions 


produced were also investigated using an analytical model. 





TIT. SUMMARY OF HEATING TECHNIQUES 


A. DIRECT CURRENT (dc) METHOD 

Localized lesions have been produced in tissue by applying a direct 
current betveen two electrodes. An active electrode is inserted into a 
predetermined location in the tissue, and an inactive electrode, or 
ground, is placed on some other part of the body. The resistivity of 
the tissue causes electrical energy to be converted into heat, thus 
increasing the local tissue temperature. When the tissue temperature 
is raised above 55°C [2], the tissue is destroyed. 

Although small and discrete lesions can be produced by the dc 
method, it has several disadvantages. Direct current follows the path 
of least resistivity and tends to have preferred paths in tissue, thus 
causing irregular lesions. Gas bubbles, a result of electrolysis, form 
around the active electrode and block the flow of current. In addition, 
dc devices suffer from polarization which can produce transients that 
stimulate muscle and nervous tissue [23]. 

Nerve stimulation does have one advantage. The probe placement 
can be checked by observing the patient's reaction to the stimulating 
effect of the current. However, other nerves, as well as those in the 
lesion, could also be stimulated. 

Herrero [14] states another disadvantage of the dc method: 

"Animals with -dc lesions usually had extensive scar tissue both around 
and infiltrated with the tissue". 

It should also be noted that the dc method depends on both the 
thermal and electrical properties of the tissue. The thermal properties 


are known within an accuracy of approximately 5% [7] and are relatively 
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constant throughout the tissue. However, the electrical properties are 
not well known, are anisotropic and tend to vary wtth location. 

One advantage of the dc method is that the electrodes are small; 
approximately two millimeters in diameter. Therefore the probes can 
be painlessly inserted and removed, and cause little damage to the 
- Surrounding tissue. 

In summary, although the dc method can produce localized lesions, 
the irregularity of the lesions, control difficulties, and possible 


Side effects of the current, greatly diminish its usefulness. 


B. RADIO-FREQUENCY (rf) CURRENT METHOD 

The radio-frequency technique produces lesions in much the same way 
as the direct current method. A long, thin, needle-like probe emits 
low power continuous wave radio-frequency current. An indifferent 
electrode is placed on some other part of the body and acts as a ground. 
As current is passed between the electrodes, the resistivity of the 
tissue causes electrical energy to be converted to tissue internal 
energy, thus increasing tne local tissue temperature. As in the dc 
method, lesions are a function of the tissue's thermal and electrical 
properties. 

Gengler [10] presents a detailed analysis of the rf technique. 

In his summary of the method he states: "The radio-frequency method 
allows controlled and predictable lesions to be produced within safe 
limits for physiological and clinical applications." 

Unlike direct current, rf current has no preferred path in tissue 
and no stimulating .effect. However, low-frequency alternating current 
does have a stimulating effect in tissue. Therefore, the probe location 
can be checked by passing low-frequency current between the electrodes 


and observing the patient's response. 
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Brodkey and others [3] discuss another way of assuring proper 
probe location. They found that a temperature rangeof 40-49°C reversibly 
blocks nervous function. Thus, after positioning the probe, and raising 
the local tissue temperature to the reversible range, the patient's 
response can indicate if the probe 1s positioned correctly. 

Like the de electrode, the radio-frequency probe is small, can be 
painlessly inserted and removed, and causes little damage to the 
Surrounding tissue. 

The rf method does have several disadvantages. The local tissue 
temperature may exceed 100°C causing boiling, gas formation, carboniza- 
tion of the tissue, and adherence of coaguluum to the probe [2]. Also, 
defects in the probe insulation have been Known to cause irregular 
lesions along the probe tract. However, these difficulties can be 
avoided by carefully controlling the various physical parameters used 
in the rf technique. 

In summary, the radio-frequency method can be used to produce well- 
circumscribed lesions and the probe location can be checked by two 
methods. However, as in the dc method, lesions are still a function of 


both the thermal and electrical properties of the tissue. 


C. HIGH RESISTANCE HEATING PROBE METHOD 

In this method, a high resistance electrode is inserted into the 
tissue. Current is passed through the electrode and the resistance of 
the electrode.causes the probe temperature to increase. Heat 1s 
conducted away from the probe into the surrounding tissue. A region 
of destroyed tissue 1s formed wnen the local tissue temperature 15s 


raised above approximately 55°C [2]. 





Theoretically, there is less uncertainty in producing lesions 
purely by heating than there is in the direct current or radio-frequency 
techniques. The dc method produces both thermal and electrolytic tissue 
destruction. As the lesion forms, it offers an increased resistance to 
the passage of current. Thus, the parameters are changing as the lesion 
forms. The rf method minimizes electrolysis but still depends on the 
electrical and thermal properties of the tissue. On the other hand, the 
high resistance probe depends only on the tissue's thermal properties. 

Carpenter and Whittier [4] attempted to produce lesions in animals 
by using a heating electrode. Although their results were poor, they 
concluded that it was due to the crude instruments employed. They 
stated: "It would seem profitable to explore this technique with a more 
finely perfected thermocautery." 

Gildenberg [11] conducted studies with a high resistance heating 
electrode and a rf probe and found that the heating probe produced less 
variable lesions than the rf probe. In addition, the rf lesions had 
Significantly more lesions along the puncture canal than lesions formed 
by the heating probe. 

As discussed earlier, a temperature range of 40-49°C can reversibly 
block nervous function. This technique can be used to check the 
position of the heating electrode. 

Since lesions formed using the resistance heating method depend 
only on the temperature field surrounding the probe, the surface 
temperature of the probe must be accurately measured and controlled. 
This has been one of the difficulties of previous resistance probe 
designs. Watkins [23] used a probe with a thermocouple placed in the 
tip. He found that the actual temperature of the probe could exceed 


the recorded temperature by several degrees if the thermocouple was not 
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positioned carefully. Thermocouples placed in the surrounding tissue 
to record the temperature can distort the temperature field and give 
erroneous readings. Gildenberg [11] oe thermocouple one 
millimeter from the heating electrode to record the tissue temperature. 
He determined that controlling the current to maintain a constant 
temperature at the thermocouple was no more accurate than maintaining 
geeonstant Current. 

Tt should also be noted, that since the lesioms produced by the 
heating probe depend only on thermal destruction of the tissue, the 
lesions are smaller than those produced by the direct current or radio- 
frequency current method. 

In summary, a properly designed high resistance heating probe can 
produce small, safe and predictable lesions. The uncertainty is less 
than the radio-frequency or direct current method since tissue destruction 


depends only on the thermal properties of the tissue. 
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TTI. HIGH RESISTANCE HEATING PROBE 


A. THEORY 
It is well known that the flow of electric current against any 
resistance is accompanied by a dissipation of electrical energy in the 


form of heat. The amount of heat generated is expressed by Joule's 


Taw: 
= eG 
QO = IV = IR (1) 
Q = heating rate 
V = voltage 
I = current 
R = resistance 


The resistance of a conductor at ordinary temperature 1s a function 
of the length and the cross-section of the conductor, varying directly 
with the length and inversely with the cross-section according to the 


following relationship: 


R= e'F (2) 
R = resistance observed 

g = length of sample 

A = cross-section area of sample 


The proportionality constant p, is called the resistivity of the 
material. Therefore, all other factors remaining equal, the material 
having the highest resistivity will generate the most heat. 

All metals exhibit some change in resistance with temperature. 
For pure metals, there is a definite and constant relationship between 
the temperature and the resistance of the metal. This characteristic 


is used in electrical-resistance thermometers where temperature 1s 


18 


indicated by measuring the change of resistance. Since electrical 
resistance measurements can be made with high accuracy on bridge 
circuits, this type of instrument can give precise temperature 
measurements. 

The theory of the resistance thermometer was applied to the design 
of the high resistance heating probe. From the above discussion, a 
given metal will have a particular resistance at a particular temperature. 
The inverse of this was used in designing the probe. That is, imposing 
a particular resistance on a metal should cause the metal to come to the 
corresponding temperature for that resistance. The resistance can 
easily be controlled by regulating the voltage and current. Thus, the 
temperature of the probe can be measured and controlled accurately 
without the use of any external measuring devices. 

The change in resistance with temperature of a metal is given by 
the following expression: 


Ry = Ry [1 a! (1 - T,)I (3) 


OO 
UW 


1 resistance at temperature T, 


AO 
tt 


9 resistance at temperature Tp. 
The constant a', is the temperature coefficient of resistivity and 
iS expressed as ohms per ohm per degree. It should be noted that, in 


general, a' iS a function of temperature. 


B. DESIGN 

Fauations (1), (2) and (3) express several factors that were con- 
sidered in designing the high resistance probe. From equations (1) 
and (2) it can be seen that resistivity should be high to generate the 
most heat. Equation (3) shows that the temperature coefficient of 


resistivity should also be high so that the change in resistance with 
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temperature will be large and therefore easier to measure and control. 


Ideally, a’ should be constant over the range of temperatures to be 
measured. 

Several metals were investigated. Their properties are listed in 
Table I. In addition, the properties of thermistors were also 
investigated. Although they are small, rugged, can be made in any shape 
and have a long life, their change in resistance with temperature is 


extremely nonlinear. Since this makes temperature control difficult, 


thermistors were not considered further. 


TABLE I 
Resistivity Coefficient 

Metal (auecuce) of Resistivity Constant 
Platinum 10 microhm/cm .0038 ohm/ohm/°C 
Nickel 7.8 .006 
Nichrome 100 .0004 
Copper hie .0038 
Aluminum Ges . 0039 
Balco Lo .0045 25-100°C 
Tron 10 .005 


After comparing the properties of the metals in Table I, Balco was 
chosen. Balco is the commercial name of an alloy, consisting of /0% 
nickel and 30% iron, that has the properties that are desirable ina 
heating electrode. Its temperature coefficient of resistivity 1s large 
(a' = 19.9 microhm/cm) and nearly constant between 25°C and 100°C [24]. 

Electrode insulation was an important consideration in the probe 
design, since a voltage leak in the probe can cause tissue electrolysis. 


Balco 1S available with a standard enamel insulation. The enamel is 
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made from an oil modified expoxy based resin and has a continuous duty 
operating temperature of 150°C [24]. 

Various shaped heating probes were considered. However, to 
facilitate the construction of the probe, a cylindrical shape was chosen. 
A probe length to diameter ratio (L/D) of ten to one was used to reduce 
end effects and to ensure good experimental data. 

In order to increase the extent of the temperature field in the 
test medium, thus increasing the accuracy of the experimental data, the 
probe was made approximately four times larger than an actual surgical 
probe. 

The heating electrode (see Figures 1 and 2) consists of a skeleton 
of 0.3894 centimeter (cm) diameter glass tubing. Two 0.038 cm copper 
wires were inserted in the tubing. One wire was passed out the tip 
and high resistance 0.0068 cm Balco wire was attached. The Balco wire 
was tightly wound back along the glass tubing for approximately 3.9 cm. 
It was then led back into the tubing through a small hole and attached 
to the second copper wire. The holes at the tip and side of the glass 
tubing were closed with expoxy. The probe was then given a thin 
coating of General Electric No. 7031 Adhesive and Insultating Varnish 
thinned with Xylene to ensure against any insulation defects in the 


Balco wire. The outside diameter of the completed probe was 0.4145 cm. 


C. TEMPERATURE CONTROL 

Precision resistance measurements can be made by using a Wheatstone 
bridge. The bridge (see Figure 3) consists of four resistance arms, a 
voltage source and a detector. When A is to B, as X is to the standard, 
there will be no potential difference bewteen points d and b. Therefore, 


whenever the detector reads zero voltage from point d to b, the bridge 
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Figure 1. Construction of Probe. 
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Figure 3. Wheatstone Bridge. 
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is in balance and the unknown resistance (X) is obtained by using the 
relationship: 
A 


=a B ) x (Standard) (4) 


The probe was calibrated using a Rosemont variable temperature oi] 
bath and commutating bridge. The calibration curve is shown in Figure 
.5. Because of the experimental apparatus, a probe temperature between 
45°C and 50°C was desired. The reason for this choice will be discussed 
Jater. It should be noted, however, that in actual use, a probe 
temperature between 85°C to 95°C would be used. 

As discussed earlier, the flow of current against a resistance 
causes heat to be generated. There is also a change in resistance with 
temperature. Since the temperature of the probe must be controlled 
precisely, it was necessary to ensure that the bridge resistances would 
not change as heat was generated in them. Therefore, 56.2 ohm resistors 
rated at five watts were used for resistances A and B. Referring to 
equation (4), it can be seen that since the ratio of A to B is one, the 
Standard resistance will also be the resistance of the unknown (probe). 

The standard resistance was wound using Evanohm wire. Evanohm 1s 
a commercial wire used in precision resistance standards and is an alloy 
of 75% nickel, 20% chromium, 2.5% aluminum and 2.5% copper. It has a 
low temperature coefficient of resistance over a wide range of 
temperatures [24]. The completed Evanohm standard had a resistance of 
356.33 ohms which corresponded to a probe temperature of 48.5°C (see 
Figure 5). 

The final probe circuit (see Figure 4) consisted of the two 56.2 
ohm resistors, the Evanohm standard and the probe. A Lambda Model LP 
413 FM regulated power supply was used to provide power to the probe. 


A Weston Model 1240 digital voltmeter was used as the detector. 
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Figure 4. Temperature Control Circuit. 
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Figure 5. Calibration Curve. 
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Since the probe had a high temperature coefficient of resistance, 
the probe resistance was constantly changing as it generated heat. 
However, by manually regulating the Lambda power supply the bridge could 


be kept balanced and a probe temperature of 48.5°C maintained. 
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IV. EXPERIMENTAL APPARATUS 


A. TEST MEDIUM 

A region of nervous tissue iS guite complex. However, to test the 
heating probe, a simplified medium that had the heat transfer character- 
istics of tissue was desired. Since the lesions produced by the heating 
probe are a result of thermal conduction, heat transfer by convection 
had to be negligible in the test medium. In addition, the medium had 
to be clear so that the region Surrounding the probe could be seen. 

Grayson [13] suggested that the water content is a major factor 
in determining the thermal conduction of animal tissue. Spells [22] 
investigated thermal conduction in both fluids and tissue and found 
there was a linear relationship between the conductivity and the water 
content of the tissue. The linear relationship was particularly close 
when the water content was above 50%. 

Cooper and Trezek [7] obtained results that indicated the thermal 
properties of human tissue can be predicted within 5% by measuring the 
water content and assuming the remaining portion of the tissue 1s 
composed of equal amounts of fat and protein. They also found that the 
percentage of water in human tissue varied between 71% and 84%. 

The above results indicate that an ideal test medium should have a 
high percentage of water. Gels would appear to be perfect. Gel 
consisting of gelatin mixed with an appropriate percentage of water 
would approximate tissue as a mixture of protein and water, and thermal 
properties close to those of human tissue could be expected. However, 

a characteristic of gels is that they melt at PeaRe cecil, 22.64. sSINee 
this would lead to convection effects, gel was not used as the test 


medium. 
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It was found, however, that agar-water solutions consisting of as 
little as 0.25% agar are clear, solidify at room temperature, and remain 
solid up to a temperature of approximately 85°C. In addition, the thermal 
properties of a 1%, or less, agar-water solution are almost identical 
with those of pure water, and within 10% of the thermal properties of 
most human tissues [7]. For these reasons, a mixture of 97.7% water and 


0.3% Difco Bacto-Agar was chosen as the experimental test medium. 


B. LIQUID CRYSTALS AND TEMPERATURE MEASUREMENT 

Under appropriate conditions of temperature and pressure, matter 
can exist ina solid state, isotropic liquid state or gas state. The 
majority of organic and inorganic compounds are comprised of a regular, 
three-dimensional array of molecules or ions in the solid state, j.e., 
gmeuystal. 

When the crystal lattice is heated, thermal agitation makes it 
impossible for the molecules to maintain a regular arrangement. At a 
given temperature, the crystal begins to melt and changes to a liquid. 
This change, known as the melting point, is usually very sharp and a 
characteristic of the compund. At still higher temperatures, the liquid 
boils and passes into the gas phase. 

Near the turn of the century it was discovered that some compounds 
had a well defined melting point but the melt was opaque. The opacity 
disappeared when the temperature was increased to a definite limit, 
given a true isotropic liquid. Although this intermediate phase was 
fluid to some extent, it also exhibited some characteristics of a 
crystalline solid. This phase was called a mesophase and compounds 


exhibiting this phenomenon became known as liquid crystals [12]. 
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There are three types of mesophase [12]. 

(i) smectic - molecules have a parallel orientation and are in layers 
with their axes at right angles to the layers. 

(i7) namatic - molecules have a parallel orientation of their axes. 
(111) cholesteric - molecules have a parallel senate ich and are in 
layers with their axes in the plane of the layer. 

The cholesteric mesophase is quite different from the smectic and 
nematic and possesses unique optical characteristics. Cholesteric 
liquid crystals having a certain structure exhibit dramatic changes in 
color with very small changes in temperature. This color-temperature 
relationship is constant and reappears everytime the liquid crystal 
passes through a particular temperature range. 

Dowden [9] discusses the structure of cholesteric liquid crystals. 
These crystals are derived from cholesterol or other sterol systems 
and have a basically planar structure. The axial direction in each 
successive layer is displaced slightly from that in the previous layer 
giving a spiral trace. The thickness of each layer is about 3A and 
the pitch of the helix is about 5500A . 

The property which 1s the basis of color behavior is circular 
dichroism. An incident beam of unpolarized light is split into two 
components having electrical vectors rotated in opposite directions. 
One component is transmitted and the other is scattered. The scattered 
light normally has a wavelength peak that occurs in the visible region 
of the spectrum. 

Temperature changes disturb the intermolecular forces and cause a 
shift in the molecular structure. This changes the wavelength of the 
light that is scattered. Thus, cholesteric liquid crystals exhibit 


dramatic changes in color with very small changes in temperature. 
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Since the colors scattered by liquid crystals are unique for specific 
temperatures, the quantitative measurement of temperature is possible 
to an accuracy of less than 0.1°C [8]. | 

Liquid crystals can be applied directly to a surface. However, a 
dark uniform layer is normally placed on the surface before the liquid 
crystals. This provides a Standard surface texture, intensifies the _ 
resulting colors and isolates the crystals from surface contamination. 

Liquid crystals have some application problems. Since they are 
viscous, applied film thicknesses can easily vary and result in non- 
reproducible thermal patterns. Direct exposure of the crystals to 
solvent vapors, oils, dirt and ultra-violet light can cause variations 
in their sensitivity. However, it 1s possible to encapsulate the liquid 
crystal in a tiny capsule (10 to 30 microns), thus shielding it from 
most of these detrimental effects [17]. 

Encapsulated crystals are available commercially in a variety of 
forms that cover a wide range of temperatures. Rochrome liquid crystal 
tapes manufactured by Hoffman-LaRoche Inc. were used. These crystals 
were coated on a mylar subtratum with an absorptive black background. 
The overall thickness of the crystals and mylar was less than 0.01 
inches. 

The Rochrome crystals are calibrated with an accuracy of 0.1°C and 
have a response time of less than one second [18]. Parker [19] found 
that the response time of the encapsulated crystals depends on the 
crystal thickness. For a thickness of 0.0035 inches, response time was 
of the order of 200 to 300 milliseconds. 

By placing the liquid crystal tape in the test medium, the tempera- 
ture field around the probe could be seen and temperatures throughout 


the field accurately measured. 
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Rochrome liquid crystal tape lot no. AS-/0-36 was used in the actual 
experiments. The color to temperature relationship of this tape was: 
29.9°C - red | 
31.5°C - green 
34.2°C - blue 
During preliminary experiments with the liquid crystal tapes it was 
found that the water in the test medium caused the encapsulated crystals 
to separate from the mylar backing. To prevent this from happening, the 
tapes were given a thin coating of Varathane and the edges were sealed 


with a rubber adhesive. 


Ccaicon CELLS 

Two test cells were used to study the three-dimensional temperature 
field surrounding the probe (see Figure 6). A circular container, five 
inches in diameter and four inches deep, was constructed of plexiglass 
and was used to visually observe the radial temperature field. A 
removable ring was fitted into the cell to support the liquid crystal 
tape (see Figure 7). The tape was attached to a thin sheet of mylar 
that was streched across the ring. The probe was inserted through the 
bottom of the cell and was held in place by a bushing that allowed the 
height of the probe to be adjusted. 

The second cell, a 4 5/8 x 5 1/2 x 6 inch rectangular container, 
also constructed of plexiglass, was used to observe the axial temperature 
field near the probe. The probe was inserted through the side of the 
container. A removable bracket was manufactured to hold the liquid 
crystal tape lengthwise along the probe (see Figure 7). Because of 
symmetry, the temperature field on only one side of the probe was 


measured. 
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To ensure that the container walls would not disturb the tempera- 
ture field around the probe, both containers were constructed so that 


their dimensions were large compared to the probe dimensions. 
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V. EXPERIMENTAL PROCEDURE 


The test cells were filled with the 0.30% agar-97./79% water test 
medium for the experimental runs. The probe was positioned in the 
circular test cell so that there was an equal length of probe above and 
below the liquid crystal tape. Thus, the temperature field around the 
midpoint of the probe was measured. The probe was covered by the agar- 
water medium so that the top of the probe was at least five probe 
diameters below the surface. 

All test runs were conducted at room temperature using Rochrome 
liquid crystal tape lot no. AS-/70-36 to measure the temperature field 
around the probe. As mentioned earlier, this tape had the following 
temperature to color relationship: 

200 0 er eG 
31.5°C - green 
34.2°C - blue 

A narrow strip of Series 9 Rochrome tape was placed in the test 
medium to measure the medium's ambient temperature. The temperature to 
color relationship of the Series 9 tape was: 

23.0°C - red 
24.0°C - green 
27.0°C - blue 

The ambient room temperature, normally about 23.5°C, was used as a 
rough check of the temperature obtained from the Series 9 tape. 

Before the containers were filled with the agar-water test medium, 
four reference marks, 4.5/ centimeters apart, were placed on the 


corners of the tape. Figure 8 is a top view of the circular test 
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cell showing the probe in place, the narrow strip of Series 9 tape, 
and the four reference marks. 

The resistance of the probe was controlled manually by adjusting 
the Lambda power supply to maintain zero voltage on the detector. 
Figure 9 is an overall view of the apparatus. 

Photographs of the temperature field surrounding the probe were 
taken using a Graflex 4 x 5 Press camera fitted with a Polaroid Land 
film holder. Polaroid 4 x 5 Land film Type 57, 3000 speed was used. 

Photographs were taken as various increments of time. The outer 
radius of the temperature field, corresponding to 29.9°C, and the 
reference marks were measured in each picture. Using the actual 
distance between the reference marks (4.57 cm) and the probe radius 


(0.2076 cm), the non-dimensional radius was calculated. 
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VI. THEORETICAL ANALYSIS 


A. ONE-DIMENSIONAL MODEL 
1. Metabolic and Blood Flow Effects Neglected 

In order to simulate the heat transfer process in the region 
Surrounding the probe, a one-dimensional heat transfer model must be 
developed. In this model it will be assumed that the tissue is an 
infinite homogenous medium with known conductivity, density and specific 
heat and that it is initially at a uniform temperature of TO: The 
probe is assumed to be a cylinder of infinite length and radius a with 
a constant temperature of Vs (see Figure 10). Thus, the problem is to 
solve for the temperature distribution in a region bounded by a circular 


cylinder. 


The general conduction equation is: 





l= ac 2 (5) 
or 

k = thermal conductivity 

0 = density 

Cc = specific heat 

a = « = thermal diffusivity 


It is now convenient to introduce the following nondimensional 


Variables: 
R= (distance) 
r 
0 
tT = “5- (time) 
ro 
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Figure 10. Mathematical Model of Probe. 





9 = a (temperature) 


Substituting these variables into equation (5) yields, 





Rp 99 
, oe. o 


The boundary conditions are: 
(a) at the probe surface (r = a, or R=1), T = tr and 6 = 1. 
(b) at large distances from the probe the temperature of the tissue 
remains constant; that 1s, as r>», T > ile or nondimensionally, 
as R>m, 6 > 0. 
The initial condition is at time t = 0, T = ie or 9 = 0 for all 
locations. 


Carslaw and Jaeger [5] give the closed form solution of equation 


z Be. C, (u, Ru) du 
a ee a 


© u [2 (u) + ¥8 (u)] 


g=1- 


C._(u, Ru) = Jo ey 


; (Ru) - ie (u) Jo (Ru) 


O 
where, Jo and iP are Bessel functions. Equation (7) was evaluated 
numerically by Jaeger for various values of t and R [15]. 
2. Metabolism and Blood Flow Effects Included 

Pennes [20] was the first to account for metabolic and blood 
flow effects in the heat equation. Pennes suggests that the equation 
take the form: 

ce 


ee. 00a 
r 


or b Cp (F 


b 
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where: 


my = blood flow rate per unit volume 
or = metabolic heat rate per unit volume 
Te = blood temperature 


In this model, the effect of blood flow is treated as a perfect 
heat exchanger. That 1s, it 1s assumed that the blood enters a unit 
volume of tissue at the arterial temperature, Th and acquires the 
tissue temperature as it exits the volume. 

If the metabolic heat generation rate 1S assumed to be constant, it 


can be combined with the blood flow term as: 





where 
a 
T =f7, + 
0 b My, Ch 


This corresponds to redefining the reference temperature. Cooper 
and Trezek [6] show that for brain tissue, i is on the order of 0.5°C 
higher than the arterial blood temperature, The 

Equation (8) may be normalized by defining a nondimensional "blood 


flow parameter”: 

mc, re 
Dai 

en Sena (blood flow term) 


Substituting this variable along with R, t and @ into equation (8) 





yields: 
i =) 
] oT een 
R em ae OT (9) 


The boundary conditions and initial condition are the same as were 
used for equation (6). 
Equation (9) was solved with the use of Laplace Transformations and 


the following closed form solution was obtained: 
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2 
K (VaR se oe (eae 
aa K, (YER ) = cof J a (10) 
K (VB) oF — [ute/u] [05 (u) + 2 (u)] 
C, (u, Ru) = J. (u) i (uR) - Jo (uR) Ne (u) 


Equation (10) was solved numerically for various values of R and t 
and for values of 8 = 0.0, 0.1, 0.5 and 1.0 and the results are tabulated 


in Tables II - XI. The reader is referred to Appendix A for details. 


B. TWO-DIMENSIONAL MODEL 

The one-dimensional model assumed that the probe was infinite in 
length. However, the probe consists of a constant temperature cylinder 
of finite length and a cylindrical insulated stem. The one-dimensional 
Solution does not take end effects into account. Thus, it was desired 
to have a two-dimensional solution that would include these effects. 

A digital computer program, by the code name of Trump [16], was 
used to determine the two-dimensional transient temperature field 
around the probe. 

Since TRUMP obtains a solution in dimensional form, a nodal netword 
(see Figure 11) had to be constructed for the geometric configuration 
of the probe. Due to symmetry, it was necessary to generate the field 
for only one half of the cross-sectional area. The dimensions of the 
network were chosen large enough so that the boundary conditions at 
R = 30 and X = 24.0 ae did not affect the temperature field during the 
transient period which was investigated. The network therefore 
Simulated an "infinite" field. 

The TRUMP solution gave values of 6 at each node point for various 
values of nondimensional time. The results are given in Tables XII - 


XVI. 
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VIT. RESULTS 


The preceding section gave details of the three models that were 
used; a one and two-dimensional model without blood flow and metabolic 
effects and a one-dimensional model including these effects. Numerical 
results for the one-dimensional models are given in Tables II-XI. 
Tables XII - XVI give the numerical results for the two-dimensional 
model. Values for 6 are given at various nondimensional locations (R) 
and times (rt). The experimental results, differences between the one 
and two-dimensional models, and the relative effects of blood flow and 


time will now be discussed. 


Pemex PERIMENTAL RESULTS 

Figures 12, 13, 14, 15 and 16 are representative pictures of the 
radial temperature field surrounding the probe taken at +t = 0.0, 1.0, 
Seepeoeo and 10.0. Figure 17, 18, 19, 20 and 21 depict the axial 
temperature field at t = 0.0, 1.0, 3.0, 5.0 and 10.0. The latter 
pictures also show end effects. 

It was convenient to use the following nondimensional variables in 
calculating the experimental results: 

be - T 


ge =——." (critical temperature) 
C Vy - ls 


where, le is the temperature corresponding to a given liquid crystal 
color. In the experimental runs, this color was red and corresponded to 
a temperature of 29.9°C. In a surgical application, if is the 


temperature at which tissue destruction begins. 


is 
Be tae (radius of lesion) 
a 


The experimental data and results are given in Table XVII. 
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The experimental results were compared with the two-dimensional 
model. Figure 22 is a typical graph of the nondimensional time for a 
critical temperature, C= 0.200. The solid line represents the two- 
dimensional solution and the dashed lines represent the error bound. 
The reader 1s referred to Appendix B for details of the error analysis. 

The experimental results fall within the error bounds except at 
Short times. This 1s due to the fact that the analytical solutions 
were obtained assuming a step input. That is, it was assumed that the 
probe came to its final temperature immediately. Because of equipment 
limitations this was impossible to duplicate in the experimental runs 
and it took the probe approximately 15 seconds to reach its final 
temperature. 

Figure 23 1s a graph of the nondimensional lesion radius versus 
nondimensional time for a critical temperature of 8. = 0.235. The 
solid line represents the two-dimensional solution. 

The axial temperature field 1S shown in Figure 24. Again, the 
solid lines represent the two-dimensional solution. Only one experi- 
mental run was conducted using this configuration. Although the 
eerinental isotherms exhibit the same shape as the analytical iso- 
therms, they are lagging in time. A great deal of difficulty was 
encountered in placing the liquid crystal sheet in a uniform plane next 
to the probe and, in fact, the edges of the sheet tended to sag after 
a short period in the agar-water solution. More work is needed on 
this particular aspect of the experiment since experimental errors are 


obviously large. 
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B. ONE-DIMENSIONAL MODEL 
1. Effect of Blood Flow on the Temperature Field 

Figure 25 1s a graph of the nondimensional temperature versus 
the radial location in the tissue for a nondimensional time equal to 
five (corresponds to a real time of approximately 150 seconds for the 
experimental probe). Four curves are presented: the dashed line 
represents no blood flow effect and the solid lines represent 8's of 
0.1, 0.5 and 1.0. Note from the graph that blood flow tends to decrease 
the temperature field, resulting in a slightly smaller lesion. To gain 
a feeling for the magnitude of this effect, a typical value of g will 
be calculated. A typical value for mM, Ch of blood flow in the brain is 
Oe cal/em?/sec/°C [21] and a typical value of the thermal conductivity 
of brain tissue is 0.0014 cal/cm/°C/sec [7]. Therefore, for brain 


2, where r. is the probe radius 


tissue, 8 iS approximately equal to 7 ros , 


expressed in centimeters. The experimental heating probe had a radius of 
approximately 0.2 cm; thus, 8 would be of the order of 0.28. However, 
an actual probe would be approximately one fourth the size of the experi- 
mental probe and g would then be much less than 0.1. Figure 26 shows 
that values for blood flow in this range decreased the temperature only 
Suga) y . 
2. Effect of Increased Time on Temperature Field 

Figure 26 is a graph of the nondimensional temperature versus 
nondimensional time. One can see that at short times, the rate of 
lesion growth is large, but as time increases, the lesion growth 
decreases and approaches a steady size. 

Figure 27 illustrates the effect of blood flow at various times. 
In this graph the temperature at a radius of 2.5 is plotted against time 


for 8 = 0.0, 0.1, 0.5 and 1.0. At short times the blood flow effect 
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seems to be negligible but as time is increased, the effect of blood 
flow becomes more noticeable. Since the probe radius enters the blood 
flow parameter aS a squared term, the blood flow effect would become 


an important parameter for a larger probe at longer times. 


C. TWO-DIMENSIONAL MODEL 

Results for the two-dimensional model are tabulated in Tables XII - 
XVI. The results are given for various positions along the probe 
length, where the X=0 position corresponds to the midpoint of the probe. 

The two-dimensional solution was used to compare the experimental 
results with the theoretical results. Figure 24 shows the isotherms 
along the probe length for various nondimensional times. By extending 
this graph it can be seen that the temperature field around a cylindrical 
probe has an elliptical shape. The dimensions of the ellipsoid can be 
changed by varying the length to diameter ratio of the probe. If this 
ratio is made small enough , the shape of the temperature field will 
approach that of a sphere. 

Figure 28 is a graph of the axial position (X) versus the radial 
‘location in the tissue for a nondimensional time equal to five and Cea 
0.300. The solid lines represent the two-dimensional solution for an 
L/D of ten, and the dashed line represents the one-dimensional! solution. 
Around the midpoint of the probe, the two solutions are in close 
agreement. However, as the distance from the mid-point increases, the 


one-dimensional solution 1S no longer accurate. 
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VIII. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

Previous experimental methods used to predict lesion size and shape 
have generally reduced to trial and error. The usual procedure has been 
to produce lesions in animals, vary the time and current, and analyze 
the results 

The theoretical predictions of the temperature field obtained from 
the two-dimensional heat transfer model of the probe fall well within 
the error bounds of the experimental results. Thus, it appears possible 
to analytically predict the size of lesions acdueed ait the constant 
temperature resistance probe. 

Experimental isotherms are reproducible and exhibit sharp delimita- 
tion, particularly at short times. .The probe itself is safe, easy to 
construct and simple. The temperature control circuit is small, 
accurate and easy to operate. 

All results are presented in nondimensional form. However, it 
Should be noted that both the experimental and theoretical models are 
Simplified and the results have been verified only over a narrow 
temperature range and only for a cylindrical probe with a length to 
diameter ratio of ten. 

The blood flow effect is found to be an important parameter in 
determining the temperature field in tissue if a large probe is used. 
For very small probes, the effect is marginal. Including the effect of 
blood tends to slightly decrease the temperature throughout the field 
and a steady state or equilibrium is reached at an earlier time than if 


blood flow is neglected. As time is increased, the effect of blood 
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flow becomes more predominant. Thus, if heat is applied for any sub- 


stantial period of time, the effect of blood flow is felt more strongly. 


B. RECOMMENDATIONS 

It 1s recommended that additional experimental work be carried out 
with the heating probe over higher temperature ranges and with various 
probe dimensions and L/D ratios. 

As mentioned previously, problems were encountered in attempting 
to place the liquid crystal tape uniformly along the axis of the probe. 
Accurate axial isotherms were therefore not obtained in the present set 
of experiments. Further investigation is needed on this aspect of the 
experiment. 

It 1s also recommended that a technique for automatically controll- 
ing the probe surface temperature be developed. Two highly desirable 
features of such a control system would be: 1) an ability to select a 
range of probe surface temperatures, and 2) an ability to pre-select 
the duration of the heating period. These two features, together with 
the analytical model, will allow a more accurate prediction of lesion 


size. 
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APPENDIX A: ONE-DIMENSIONAL MODEL SOLUTION 


Statement of Problem: 


A cylindrical probe of radius r 


0 


1s embedded in tissue which is 


initially at a uniform temperature ee Within the tissue there is a 


volumetric heat generation due to metabolism of strength Sin" In 


addition, there is blood flow in the medium with a flowrate Ms» 


Specific heat Ch, and temperature T,- Assume that the length to dia- 


meter ratio of the probe is large enough so that end effects may be 


neglected and assume the tissue is homogeneous and infinite in extent. 


The probe is at a constant temperature T_. 


p 


Also, assume that at an 


infinite distance from the probe, no temperature change occurs. Solve 


for the tissue temperature assuming constant properties and 


one dimensional radial heat flow. 


Notation: 
Symbol 

a 

C 


Ch 
k 


Description 


Tissue thermal diffusivity 
Tissue specific heat 

Blood specific heat 

Thermal conductivity of tissue 
Volumetric blood flowrate 
Radial location in tissue 
Probe radius 

Metabolic heat generation rate 
Time 


Tissue temperature 
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Typical Units 
cm/sec 
cal/gm/°C 
cal/gm/°C 
cal/cm/sec/°C 
gm/cm?/sec 

cm 


cm 


AL enc See 


Sec 


= 








Symbo ] Description Typical Units 


Th Arterial temperature 7 a 
VA Initial tissue temperature a 
Ty Probe temperature SC 


Nondimensional Groups: 
T-T 





—_ aT (temperature) 
Pp 0 
R = — (distance) 
= 
0 
a - (time) 
ic 
2 
fii wet 
oe (blood flow) 
k 
° . ° ° 
Me Vi m, Cp (initial temperature) 


The general heat conduction equation in cylindrical coordinates is: 


ay eo 
kK ( ar 
ak 


, 
+m, cy (I, -T) +S, = pc = (A.1) 





Substituting the above nondimensional quantities into equation 


vem) yields: 


aR) 
aR _ 20 
ta | ae a 


E 





Equation (A.2) must be solved in order to describe the temperature 
of the tissue as a function of time. In order to solve the problem, 
boundary conditions and initial conditions must be stated. They are 


as follows: 


7| 





(a) at r = ro? (at the probe surface), T = T or, in nondimensional 
format R=1,60=1, 

(b) as r + », the medium remains unaffected as r >, T > iP Orn 
in nondimensional form as R > », 6 > 0, 

ke) at t =o) = ie or, in nondimensional form at t = 0, 6 = 0. 

By using Laplace Transformations, a closed form solution to 


equation (A.2) can be obtained: 





2 
KER 5 a 8 “i C. (u, Ru) 
Oa = = ino ee Ce du (A.3) 
KO VB O (v i als oe "o (u) J 
Cy = Jo (u) Ne (uR) - Jo (uR) A (u) 


If the blood flow term is neglected, equation (A.3) reduces to: 


2 
om Tu C, (u, Ru) du 
7 ——— (A.4) 
O 


5 = = 
u [95 (u) + ¥6 (ud 


a[n 


Equation (A.4) is the closed form solution given by Carslaw and 


Jaeger [5]. This integral has been solved by Jaeger [15]. 


Solution: 
A finite increment technique was used to solve for the temperature 


field. Equation (A.3) was rewritten in the following form: 


pe 2 faye (A.5) 
0 


2 
eo) (Carer 
_ O 
ee a ay ea 
(u + G)L9§ (u) + YS (u)I 
Assume an approximation to the integral can be derived by passing 


an interpolating function through several adjacent points and then 


We 





integrating the interpolating function. Simpson's Rule uses the 


parabola as an interpolatory function between U._1> U. and u.,, where 
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the u. are evenly spaced values of the interval (u_, uA) with a 


0 


Subinterval length h = u u.. Integrating the parabolic 
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interpolating function gives: 
Sit ' 
meee! 
Jd au) du= 2 Oy, +45 # Yigg J (A.6) 
where Y. = q(u,) 


Expanding equation (A.6) over the entire interval yields, 


+ Y |] CN) 


n 
me 
J a(u) du = 5 nips + AY, + 2Y,, ee AY ny 


Computer programs 1, 2 and 3 were used to solve this problem for 
B = 0.0. Programs 4, 5, and 6 were used to solve for 68 = 0.1. To 
solve for other values of B only one card must be changed. 

Because q(u) is singular at u = 0, the lower limit of the integral 
was taken as 0.001. The error caused by this was assumed to be negli- 
gible. By varying the upper limit of the integral between 1.0 and 100.0, 
it was found that the solution coverged at approximately 8.0. Therefore, 
an upper limit of 10.0 was used. 

The Y Bessel functions involving the lower limit of the integral, | 
0.001, were solved using a separate program and calling the BESY 
Subroutine. This was done to avoid underflow and reduce computer time. 

The included programs solve the temperature field for times 
between 0.1 and 1.0 only. To change the time interval only one card 
must be changed. 

Jaeger's results [15] were used as a check on the accuracy of the 
results obtained from Programs 1, 2 and 3 when g = 0.0. The greatest 


difference between the solutions was 0.001. 


Us 





APPENDIX B: ERROR ANALYSIS OF EXPERIMENTAL RESULTS 


Three nondimensional quantities were calculated in analyzing the 


experimental results. 


ae 
@ = =——— ~~ (critical temperature) 
c.6UOUWTNd edi 
p 0 
re 
R. as (critical radius) 
0 
tT = ot (time) 
rc 
0 


Ambrosius [1] gives the following equations for the resultant 


errors (R) in computations 
(a) addition 
(Asta) + (B + 


R= + (4° + b°) 


b) 
IL 
2 


(b) subtraction 


(c) multiplication 


(A +a) (B + b) 


involving measurements: 


1 


—— 


2 = + [(Ab)* + (Ba)*]2 (B.3) 


(d) division 


(A)¢ 


1 
2 2 
(Ba/A)” + (b)" 42 (B.4) 


where A and B are two measurements and an error of + ais made in the 


measurement of A and an error of + b 1s made in the measurement of B 
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The following errors were estimated for the variables in 60: 


T =+0.1°C 
E 

i arr eae 

T. = + 0.2°C 
p 


An examination of 6. along with equation (B.4) indicates that the 
largest error in 6. occurs when i 1s the greatest, or when the quantity 
i - un is the smallest. The highest value of 7 was in test run 10. 
The resultant errors will now be estimated using data from run 10. 


ee hco.c, Clemons 0 


a Secon © eae le 
t, = 46.5 Cleedece 
Ve - le = 47°C 
V, - Un = leec 6 


Substituting these values into equation (B.2), the following errors 
are obtained: 
T. - T. = + 0.32°C 
T. - Ue + oe 36 ae 


From equation (B.4) the resultant error in O. Se 
R= + 0.025 
or, the percentage error in he 12.5%. This error would be reduced if 
a higher probe temperature (T,) or a higher critical temperature ry 
were used. 

An analysis of the error in computing Re and t is more difficult 
Since the variables in these quantities change with time. It was 
estimated that at short times the error in bothRcand +t could be of the 
order of 10%. However, this error decreases significantly as time is 


Increased. 
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